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Abstract 

Using the mass balance equations for chemical reactions, we show how the system relaxes towards 
a steady state in and out of the Onsager region. In the chemical affinities space, after fast transients, 
the relaxation process is a straight line when operating in the Onsager region, while out of this 
regime, the evolution of the system is such that the projections of the evolution equations for the 
forces and the shortest path on the flows coincide. 

For spatially-extended systems, similar results are valid for the evolution of the thermodynamic 
mode (i.e., the mode with wave- number k = 0). 

These results allow us to obtain the expression for the afhne connection of the space covered by 
the thermodynamic forces, close to the steady states. Through the afhne connection, the nonlinear 
closure equations are derived. 
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I. INTRODUCTION 



It is known that the evolution equations for the moments of the distribution functions 
we are confronted with the familiar difficulty appearing so often in statistical physics: an 
infinity hierarchy of moment equations. When we truncate the hierarchy, the resulting set 
of dynamical moments equations are not closed. Thus, a theory for the closure equations, 
(or the transport theory) should be formulated and combined with the truncated moment 
equations. A thermodynamical field theory (TFT) has been established addressing the 



closure problem for thermodynamic systems out of Onsager's region [1]. This theory allows 
to determine the nonlinear corrections to the Onsager theory so that the thermodynamic 
theorems for systems far from equilibrium are respected. Clearly, transport theory without a 
knowledge of microscopic dynamical laws can not be developed. Transport theory is only but 
an aspect of non-equilibrium statistical mechanics, which provides the link between micro 
and macro- levels of descriptions. This link appears indirectly in the "unperturbed" Onsager 
coefficients used as an input in the equations. The Onsager coefficients, which depend on 
the material under consideration, have to be calculated in the usual way by kinetic theory. 
The nonlinear contributions are derived starting from the De Donder-Prigogine principle 
[also refereed to as the Thermodynamic Covariant Principle (TCP)] stating that: "Thermo- 
dynamic systems are thermodynamically equivalent if, under transformation of fluxes and 
forces, the bilinear form of the entropy production remains unaltered". Flux-forces trans- 
formations leaving invariant the expression of the entropy production, are referred to as the 
Thermodynamic Coordinate Transformations (TCT). The TCP provides the possibility of 
defining an abstract space, the thermodynamical forces space, covered by the r independent 
thermodynamic forces X z , whose metric is identified with the symmetric part of the trans- 
port matrix. The law of evolution is not the dynamical law of particle motion, or the set 
i.e., of two-fluid macroscopic equations of plasma dynamics. It is rather the evolution of the 
thermodynamic forces of the system in the thermodynamical forces space. 
The nonlinear transport equations have been successively derived from two assumptions: 

1. The thermodynamic theorems valid when a generic thermodynamic system relaxes out 
of equilibrium are satisfied. In particular, the Glansdorff-Prigogine Universal Criterion 
of Evolution (UCE)] is fulfilled; 
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2. There exists a thermodynamic action, scalar under TCT, which is stationary for general 
variations in the transport coefficients and the affine connection of the thermodynam- 
ical forces space. 

However, the determination of the thermodynamic action requires the knowledge of the 
affine connection. As we shall see, this expression can be deduced by analyzing several 
examples of relaxation. 

This theory has been successfully applied for studying transport processes in non equilibrium 
systems such as magnetically confined plasmas, materials submitted to temperature and 
electric potential gradients or chemical reactions. In particular, in tokamak-plasmas, the 
thermodynamic forces and the conjugate flows are the generalized frictions and the Hermitian 
moments, respectively. Hence, we can determine the particle fluxes (electrons and ions) 
and energy losses as well as the (nonlinear) particle distribution functions. This has been 

n n 

accomplished in Ref. |2| and in the paper reported in the footnote [15| . 

In this paper, starting from the mass balance equations applied to chemical reaction-diffusion 

systems, we establish the following analytical results: 

(i) After fast transients, in the Onsager region of the thermodynamical forces space, a 

homogeneous system relaxes towards a stable steady state along a straight line. 

(ii) After fast transients, out of the Onsager region, a homogeneous system where r chem- 
ical reactions take place simultaneously, relaxes towards a stable steady state such 
that Yli=i J% M l (X, g) = YJi=i Ji M % (X s .p., o) where U l (X, g) and U l (X s p , g) are the 
evolution equation for the chemical affinities and the shortest path in the chemical 
affinities space, respectively. Jj denotes the vector of the chemical reaction rates. The 
trajectory traced out by the system and the shortest past are parametrized by g. This 
parameter is defined in the Subsection IIII Bl [see Eq. (155])]. 

Similar theorems for the relaxation of the thermodynamic mode (i.e., the mode with wave- 
number k = 0) to the steady-state, can also be derived for spatially-extended systems. 
These theorems allow us to obtain the expression for the affine connection of the space cov- 
ered by the thermodynamic forces, close to the steady-states. Through the affine connection, 
the nonlinear transport relations are derived. 

The paper is organized as follows. In Section [Til t ne simplest cases of isothermal and 
homogeneous chemical systems, relaxing in the Onsager region, are analyzed. A detailed 



study of chemical systems out of the linear region can be found in Section IIHI Since we 
are interested in the solution close to the stationary states, all expressions and matrices, 
entering in the calculations, are be developed in terms of a small parameter e measuring the 
"distance" of the system from the steady state. We conclude in Section fllVp . 
Such examples enable us to determine the expression of the affine connection for the thermo- 
dynamical forces space, close the non-equilibrium steady-states. In Ref. [1], the nonlinear 
closure equations have been derived from assumption 2. 



II. CHEMICAL REACTIONS IN THE ONSAGER REGIME 



This section is devoted to the analysis of chemical systems relaxing towards a steady- 
state in Onsager's region of the thermodynamical forces space (i.e., the space of the chemical 
affinities) jjij]. For the sake of simplicity, we start by considering the relaxation of isothermal 
and homogeneous chemical reactions. We begin by studying in detail the relaxation of a 
triangular chemical system, maintained out of equilibrium by keeping constant the concen- 
trations of two (among the three) chemical species. Our analysis is successively generalized 
by considering systems where r independent chemical reactions take place simultaneously. 
Spatially-extended thermodynamic systems will be analyzed in the next section. 



A. Unimolecular Triangular Chemical Reactions 



Let us consider an ideal homogeneous system where the chemical reactions involve three 
substances B, C and D obeying to the following chemical scheme 4| 

B ^ C 
C^D 
D^B 

Indicating with A iy Ji and Ci the chemical affinities, the chemical reaction rate of reaction 
i, and the concentrations of species i, respectively, we have 



Ax 
A 3 



Hb - He 
He — Hd 
Hd — Hb 



■h 

J 3 



t 1 C B 
t 2 C c 



(1) 







~tiC B 


= t 1 C c 








= t 2 c D 











where ki are the direct (— >) and the inverse («— ) kinetic constants of the reactions and fii 
the chemical potentials. Note that, from Eqs ([T]) we find A\ + A 2 + A 3 = 0, so only two 
affinities, say Ai and A 2 , are independent variables. At equilibrium, each chemical reaction 
rate vanishes, so the equilibrium concentrations Cb,c,d satisfy the conditions 



(2) 



We maintain the system out of equilibrium by keeping constant the concentrations of species 
B and C: 

Cb = Cb = const. 
C c = C c (l + 5) = const. ; (5 < V 
The dimensionless parameter 5 measures the distance from equilibrium. This parameter is 
supposed to be small so that the linear analysis applies. Due to constraints Eqs ([3]), the 
concentration of species D varies. Its variation can be described by introducing a small 
dimensionless parameter e 

C D = C D (l+e) (4) 

The mass-balance equation for C B reads 
dC D 



(3) 



dt 

or, in terms of parameter e 



J2 



ds 
dt 



[k 2 + ~£ 3 )e + *k 2 5 



This system admits the non equilibrium steady-state given by 



£st. 



t 2 

— =7*"' 

k 2 + k 3 



(5) 



(6) 



(?) 



Close to equilibrium, we can write 





k yCc A x 
R T 




A L = 


< J 2 - 


k 2 C C A 2 
R T 


and < 


A 2 = 


J3- 

V 


k iC D A 3 
R T 







RT(S 
RTe 



(8) 



-{Ai + A 2 ) = -RT5 

where R is the gas constant and T is temperature of the system. By performing the linear 
transformations 

( K = A % 
J[ = J l -h (* = 1,2) 



(9) 



we obtain 

2 

ji=j2 L * xj ( 10 ) 

3=1 

2 3 

l E j m = h E ( or > ^ = * ) 



1=1 " 



where a indicates the entropy production. The expressions for the Onsager matrix Ly and 
the thermodynamic forces X* are 

1 (txCc + tsCo ~£ 3 C D \ jx 1 = ^ = R(5-e) 

^ \ ^ 3 C D k 2 C c +k 3 C D J [X 2 = ^ = Re 

Since the entropy production remains unaltered under transformations ([9]), according to the 
general formulation, the systems (Aj, Jj) and (AJ, J-) are thermodynamically equivalent jlj], 
5], We can check that the thermodynamic forces satisfy the evolution equations 
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r// " 



dX % 

J2 A j Xj + bi (i = 1, 2) with 

3=1 



A) = -(£ + ^ J j ; 6* = + t) ^ ^ j (12) 

The relaxation of the thermodynamic forces in terms of the arc-parameter q is obtained. 
From the definition of <j 

2 

dq 2 = L ij dX i dX j (13) 



we find 



where 



X x = R{8-e si )-L- x l\<;-l) 
X 2 = Re st + L-y 2 (<;-l) 



(14) 



L = {k2+ R kz) C D ; I 2 = J2 LiAXl - Xi t ){Xl - Xj.) (15) 

ij=l 

Xq are the initial conditions and = — £ s t)i R^st.]- As seen from Eq. ( IT4"1) . in the 
Onsager region of the thermodynamic space, the system relaxes to the steady-state along a 
straight line. 

Let us now quote the Glansdorff- Prigogine dissipative quantity V . Noting that 

L l/2 

e-t -e = — (Z- s ) (16) 
6 



after a little algebra, we have 



i=l 
2 



-^L- 1 / 2 ^. -£) = -(/ -?)<0 



v = E J ^ = p !l = -( ^ 2 + ^ s)(/ - ' )2 K ~ (/ " " )2 - (17) 

n 

in agreement with the general results reported in Ref. [1]. It is worthwhile to mention that 
the compact formalism provided by the differential geometry allows to obtain results (|Tj 
directly, in only one line of calculation. 



B. Several Independent Chemical Reactions 

Let us now analyze the relaxation of r coupled chemical reactions taking place simultaneously 
in a uniform system. In the Onsager region, the balance equations for the n independent 
thermodynamic forces X a read [see Eqs (IT2|) ] 

dX a r 

- W = Y. A T> XP + ba (<*,/3 = l,...,r) (18) 

a/3 

where is a matrix independent of the thermodynamic forces (but it may depend on time) 
and b a is a constant vector. It is useful to write the solutions of Eqs (}I8"]) as 

r 

X Q (t) = - V A a Bst + f a (t) with lim f Q (t) = (19) 

' * r ■ t— S> + 00 

/8=1 

where A Q st is the inverse of matrix A^j evaluated when t — >■ +oo and the steady-states, , 
satisfy the condition 

r r 

= - E with ^, = E ^a,,^, = (20) 

,3=1 A=l 

where 5^ denoting the Kronecker delta. The expression of the arc-parameter reads 



1 r+oo 1 

d? = L a pdX a dX? V L af3 (t)f<*(t)fP(t) dt = F(t) (21) 

0,(8=1 * \ 0,(9=1 

where " dot" indicates the time derivative. The second principle of thermodynamics ensures 
that dq 2 > 0. Substituting Eq. ( 12T1) into Eqs ( fT9j) . close to the steady-state, e = <C 1, 

we have 

r 

= " H A %t P + &?(? - + - O 2 + • • • (&? = const.) (22) 



Now, taking into account the symmetry 



(?-*)-► -(?-0 



(X Q -X°)^-(X Q -X«) 



(23) 



we finally obtain 



r 



= - ^ + 6?( ? - /) + 0(e 3 ) 



(24) 



/3=i 



We may then summarize our results as follows: In the Onsager region of the thermodynamical 
forces space, after fast transients, the thermodynamic system relaxes towards a stable steady 
state along a straight line. 

ILL RELAXATION IN THE NONLINEAR REGION 

In the following sections we shall be concerned with the relaxation of chemical reactions 
towards a steady-state in the nonlinear region of the thermodynamical forces space i.e., 
out of Onsager's region. We first study the relaxation of isothermal and homogeneous 
chemical reactions and then we shall extend our analysis to spatially-distributed chemical 
systems. The study of these examples will enable us to establish the validity of the statements 
mentioned in the introduction. 

A. Chemical Reactions in the Nonlinear Region 

In the first part of this section, we shall prove that, close to the steady-state, after 
fast transients, the Glansdorff-Prigogine dissipative quantity, computed through the bal- 
ance equations coincides with the one related to the system relaxing along the shortest path 
in the chemical affinities space. In the second part, we shall show the validity of the state- 
ment (ii) formulated in the Introduction. After fast transients, near the steady state, the 
projections of the evolution equation for the thermodynamic forces and the shortest path 
on the thermodynamic flows on the thermodynamic flows coincide. 

Let us consider a system consisting of n components amongst which r chemical reactions 
are possible. The (local) rate of change of the mass of component a is 




(25) 



i=l 
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where p a is the density (mass per unit volume) of a, v Q is the velocity of a, and v % a Ji the 
production of a per unit volume in the i th chemical reaction. The quantity v % a divided by 
the molecular mass M a of component k is proportional to the stoichiometric coefficient with 
which a appears in the chemical reaction i. The coefficients z/ are counted positive when 
components a appear in the second, negative when they appear in the first member of the 
reaction equations. Since the mass law is conserved in each separate chemical reaction we 
have 

n 

£X = (i = l,2,...,r) (26) 

a=l 

The closure equations for Eqs (125|) read 

r 

J t = J2nM)A 3 (* = l l 2,...,r) (27) 
i=i 

where we have introduced the chemical affinities A4 of the reaction % [i = 1,2, •• • , r) and 
the transport coefficients (which may depend on the chemical affinities). In general, ry 
is an asymmetric matrix. The chemical affinities are linked to the chemical potential of 
component a, denoted by /i fe , through the De Donder relation [3] 

n 1 n 

A = or X* = --J2»> a (< = l,2,...,r) (28) 

a=l 0=1 

where X 1 denote the thermodynamic force of the reaction i (X 1 = Ai/T). We shall now 
analyze the relaxation towards stationary states. We shall firstly study the simpler case 
of isothermal and homogeneous chemical systems and, successively, we shall extend our 
treatment to space-dependent chemical systems. 

B. Isothermal and Homogeneous Chemical Reactions 

In this case the mass balance equations, Eqs (125|) . simplify to 

r 

Pa = J2< Ji (a = l,2,--- ,n) (29) 
i=i 

From Eq. fl27H29p . we obtain the equation for the thermodynamic forces 

r r 

X 1 = - M ji {X)Jj = -J2 M^(X)r j7 (X)X' where (30) 
i=i j,i=i 

M*{X)= E^K^-(^) (M = l,2,.--,r) 

a,B=l 



By inverting matrix M l i(X) (see also the footnote [16|]), the equations for the thermodynamic 
forces can be brought into the form 

r r 

J i = ~Y^ Mij {X)X j = nj {X)X j (i = 1, 2, • • • , r) (out of the steady state) (31) 

3=1 3=1 

EI=i <HXst) = E- i= i KM l3 {X s t)X% t . = (a = 1, 2, • • • , n) 

(at the steady state) 

Z{ J=1 M tJ (X st .)X% t X^\ st . = 

where M{j(X) = My denotes the inverse of the matrix M lj i.e., Y7i=i M u Mji = 5*. The 
Glansdorff- Prigogine dissipative quantity V is immediately obtained from Eq. ( 13TT) 



V = J ^ X ' 1 = - Yl Mij( x ) xixj < ( 32 ) 

Being My(X) a positive-definite matrix (for the demonstration see, for example, Ref. js]), 
the dissipative quantity V is always negative throughout the evolution of the system. The 
inequality is only saturated at the steady-state. In literature, inequality Eq. (|32|) is referred 
to as the Universal Criterion of Evolution (UCE) Jq]. 

In Ref. [l| it is demonstrated the validity of the following theorem: // a thermodynamic 
system relaxes towards a steady-state along the shortest path in the space of the thermody- 
namical forces then the UCE is automatically satisfied and 



where q denotes the arc-parameter, defined as 

r 

dq 2 = T ij dX i dX j > (34) 

i,j=l 

and g is a convenient parameter allowing to cast the equation for the shortest path in the 
space of the thermodynamic forces into the form 

d 2 X i . dX j dX l 

= » (. = l,2,-,r) (35) 

3,1=1 

Parameter g increases monotonically as the thermodynamic system evolves in time. It can 
also be chosen so that it vanishes when the thermodynamic system begins to evolve and it 
takes the (positive) value, say I, when the system reaches the steady-state denotes 
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the affine connection of the chemical affinities space jl|. For easy reference, the expression 
of T* 7 is reported in the Appendix of this paper. We want to check that, after fast transients, 
close to the steady state, Eq. (1521) reduces indeed to the general expression, Eq. (155|) . To 
this purpose, we re-consider the forces equation, Eq. (1501) . the closure equation, Eq. (1271) . 
and constraint (1261) . Defining Y l = X 1 — X l st , close to the steady-state X st _, we may write 



M ji (X) = M ji (Y) + M3\X st ) ■ Ji{X) = Ji{Y) + Ji{X st ) 



(36) 



Taking into account that 



3=1 3=1 



in terms of the variable Y, we obtain 



Y< = - Y.]=iW ji {Y) + M^X^ y.^Y) (i = 1, 2, • • • , r) 
y*(+oo) = (i.e.,X*(+oo) = X^) 



(i = 1,2,-.. ,r) 



a=l 



We note that Eqs (1381) are manifestly invariant under transformations 



t -> -t 

yi _^ _yi 



or 



yi _^ _y< 



^1 — >■ — f„ 



with 



(37) 



(38) 



(39) 



Recalling the definition of matrix M u [see Eq. (I30p ]. we also find that, under Eqs 
variables Jj, fy and M u transform as 



(g-l)-> -(e-i) 

yi _^ _yi 

v\, — > —v\ 



< 



MY) -)■ -Ji(y) 
M ii (y) + M, i (x st ) 



(40) 



— )• 



M ii (y) + iw ii (x st ) 

Hence, close to the steady-state [i.e., for e = (g — l)/l <C 1], we may adopt the following 
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expansions (see also the footnote [1 



F* = X* - Xi t = XfaiQ -l) + X} 3) (g - If + C>(e 



(3)< 



^ = f ! = ^i)+3^ 3) (^-0 2 + O( e 4 ) 



^00 = % PO = r v( Y + x st) = Tij{Xst.) + rg\g - If + 0(e 4 ) 
MY) = J t {X) - UX st ) = jj% -l) + 4 3) (g - If + 0(e 5 ) 

MV(Y) + MV(X st ) = AfV(X) = M^_ x) {q - I)- 1 + M^(g - I) + M^(g - If + 0(e 5 ) 

( Ma = M§\q -l)+ M\?{Q - Z) 3 + 0(e 5 ) 

(41) 

By inserting expansions Eqs (jUJ) into the mass balance equation, Eqs fl38l) . at the first order 
we obtain 

r 



or, by inverting the matrix M/ 



(-i) 



(i) 



^r,,(X s ,)^ 1} + 0(6 2 



(43) 



Multiplying both sides of Eq. (j43H by (^ — Z)^Qi) an d summing over the index "i", we find 

r 

= f ^r^X^X^X^dg + O^) 



(1) 



(44) 



However, from expansions Eqs (EH), we have 



^■(X st .) = % - + 0(e 2 ) 



(45) 



M ) ( e -0 = ^i + O(e 3 



Eq. (144 p may then be rewritten as 



" V M ^ 
® 2—* lj dg dg 



i , r 



E 



dX i dXi 
' dg dg 



dg + 0(e 3 ) = J 1 [j^fdg + 0(e 3 ) (46) 



where identity Eq. ( 134")) has been taken into account. Close to the steady-state, we finally 
obtain 
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in agreement with Eq. ( 133|) and Ref. [lj]. 

We are now interested in the evolution of the system near the steady-state. To this end, 
we reconsider the equation for the thermodynamic forces, Eq. ( 131) I) , rewritten in terms of 
parameter g 

dX i r 

g =-J2M^(X)J J (< = l,2,... | r) (48) 
dg 3=1 

We note that this balance equation transforms in a covariant way under the following ther- 
modynamic coordinate transformations (TCT) [lj] 

x h = y dX^ xj 
Indeed, by inserting transformations ( 149]) into the Eq. (148]) we find 

" i=l 2,9=1 J=l 

Hence, matrix M' 3 transforms like a thermodynamic tensor of second rank. The TCT 
are the most general forces-transformations leaving invariant the expression of the entropy 
production and the Glansdorff-Prigogine dissipative quantity V [lj]. V can be obtained by 
deriving both sides of Eq. ( I48p with respect to parameter g. However, this operation should 
be performed in such a way to preserve the covariance under TCT. For this, we undertake the 
thermodynamic covariant differentiation along a curve of both sides of Eq. (|48p [L], getting 

rft Y i _ r _ A Y3 A Y l _ r _ _ r _ 

^ + E r ^^ = - £ (e- l M«4),, j-'M'%+r' £ r5,M*M*v, (5i) 

U 3,1=1 U ^ j,l,p=l j,l,P,q=l 

where comma (, ) stands for partial differentiation with respect to the thermodynamic forces. 

Eq. (15T]) may also be rewritten as 

rP X i r rlXi dX l 
U\X, g) = — + J2 V ^ — — + ki{X) = U ^ X ^ $ + = where (52) 

3,1=1 

i ,i / \ d 2 X l s. — ■ dX^ dX l , 

"'(^^^^r+j:^, — — and 

r r 

h\x)= J2(e~ lMij Jj),ie~ lMlp J P -e~ 2 ^W p m^j p j q 

j,i,p=i j,i,p,q=i 
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By using expansions ( )4T|) . close to the steady-state, we get (see Appendix) 

r 

h i (X)r^O(e) and J] Jih\X) ~ 0(e 2 ) (53) 



i=l 



Near the steady-states, Eq. (|52|) satisfies the UCE. Indeed, multiplying Eq. (|52|) with the 
thermodynamic flows Jj and summing, we obtain 

: dX' l dX l 

- \ T \ r , X • 

dg 



JJ do do 



do do 



dX l dX l 
dg dg 



0(e 2 ) = 



(54) 



where P = YH=i Ji^f~ anc ^ Eq. (E3J) has been taken into account. In the Appendix, it is 
shown the validity of the following identity 



x ^ T ^ dX i dX l _ ^ v - dX l dX l 
^ Jj il ~dJ~dJ ~ ^ TiU 1g~^g~ 

i,j,l=l i,j,l=l 



(55) 



From this equation and Eq. (1341) . we obtain the expression 





dg \dg< 

Integrating from the initial condition to the steady-state, we get 



1 ' dc \ 2 

dg + O(e 3 )>0 



(56) 



dg 



However, gP(X l st ) = 0. Indeed 

r r 



(57) 



i=l 



71 r v 

E[(E"wM(E^f 



9 ///* 







(58) 



a,/3=l i=l j=l 

where Eq. (13TT) has been taken into account. So, after fast transients, we finally obtain 



i=l 



dA i ~ /7tfc\ 2 , 



(59) 



The inequality is only saturated at the steady-state. The main conclusion of our analysis is 
thus: After fast transients, in the chemical affinities space, the system relaxes towards the 



14 



steady- state so that 

r r 

JiU\X, q) = J2 JM\Xs. V am g) + 0(e 2 ) = or 

1=1 8=1 
r r 

JiU\X, g)-J2 ■hU\X s . vath , g) ~ 0(e 2 ) (60) 

i=l i=l 

Moreover, the Glansdorff-Prigogine dissipative quantity V, differs from the one evaluated 
assuming that the system relaxes along the shortest path, by a factor of the order 0(e 3 ). We 
have also verified that, during the relaxation, the UCE is satisfied. Note that in Onsager's 
region, g oc <; (see Ref. and dr^/ds = 0. Hence, Eq. (I5T1) reduces to 

<kX + O(e) or X { =Xi L + b\^-l)+0(e 3 ) (61) 



d? 

In agreement with Eq. (|2 



C. Spatially-extended chemical systems 

Previous results can easily be generalized for spacial chemical systems. In this section we 
shall show a general result valid for the relaxation of the thermodynamical mode (i.e., the 
mode with wave-number k = 0). For this we shall adopt the following strategy: firstly 
we develop the space-time dependent thermodynamic forces (X 1 ), flows (J7i) and transport 
coefficients in (spatial) Fourier's series. Then, we perform the same calculations as in the 
homogeneous case by t ak ,n g into account the P nnc ¥e Q (see the footnote Q). 



Close to the steady-state, we have [l] 

dXi >A - dX j 

; — — dv ~ V > 

j=l * j=l 



J{k)j(t) = — J Jj(r,t)exp(-ik-r)dv 

XLJt) = l~ r I X j (r, t) exp(-ik • r)dv (62) 



dv denotes a (spatial) volume element of the system, and the integration is over the entire 
space V occupied by the system in question. Suffix "(0)" denotes the thermodynamic mode. 
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The evolution equations for p(k) a reads 

r 

P(k) a = - k - (P(k)a*v (k )) + J2 u ^y (a = l,2,...,n) (63) 

3=1 

n 

^(k) = -XX^(k) (i = l,2,...,r) 

ct=l 

where symbol " *" indicates the convolution and ft a = fi a /T. At the thermodynamic mode 
we have 

r 

P(P)° = ^2^(0)3 + O.T.(k,k') k ^kYo (« = l,2,...,n) (64) 

3=1 

n 

^(o) = ~ ^(o) (* = l,2,...,r) 

a=l 

where O.T.(k, k')k^kVo stands for oi/ier contributions from different wave-numbers. Due to 
the slaving principle, near the steady-state, the evolution equations for f>(o)a reduce to 

r 

P(o)a = ^2^aJ(o)j (a = l,2,...,n) (65) 

3=1 

n 

^(o) = -XX%)) (« = l,2,...,r) 

By performing now the same calculations as for the homogeneous case, we arrive to the 
following final result: the thermodynamic mode (k = 0) relaxes to the steady-state according 
to the law 

r r 

J( 0)j U j (X i0) , q) = J2 J{o)jU j {X {0)s . path , q) = + 0(e 2 ) (66) 

3=1 3=1 



IV. CONCLUSIONS 



The Universal Criterion of Evolution theorem has been demonstrated by Glansdorff and 
Prigogine through the balance equations. From this theorem and by analyzing examples 
of relaxation of chemical systems close to steady-states, we may obtain the expression for 
the affine connection of the thermodynamical forces space. The construction of the correct 
affine connection is made "step by step". At the first stage, we suppose for simplicity that 
the matrix of the transport coefficients is symmetric i.e., Ty = gfy. We multiply Eq. ( I48p by 
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gi P and sum over the index p. By taking the derivative, with respect to parameter g, of both 
sides of the resulting equation, we find 

_ r _ J2 yl 1 _ r _ rlY l rl YP _ r _ / \ 

E^4E(^^.,')^^ = E (i->«M*J,)i->M"J, (67) 



dp 2 2 ^— ' ' dp dp 

1=1 * l,p=l * * j,l,p,q=l 



From Eq. (IA13I) (see the Appendix) we have 

r r 

(g-'guM^J,) jf x M^j q = -6 afx^ie -i) + o(e 3 ) (68) 

j,l,p,q=l - j=l 

So, close to the steady-state, we obtain 

Ed?X^ 1 ^— \ / v dX^ dX p . . f„n\ 

ffii ^ + 2^ (SU, P + 9i P ,i) -^-^ = O(e) (69) 
i=i K i, P =i K K 

At the lowest order, previous equation may be rewritten as 

d 2 X L 1 <A iU ,dX l dX p . 

¥ + ^ sJ ^ + ^%17 = ° (e) (70) 
j,i, P =i 

where g tp denotes the inverse of the matrix gj p (i.e., Z^J v= \g %v gjp = ty). It is not diffi- 
cult to show that Eq. (ITDl) satisfies the UCE but it is not covariant under the TCT. This 
because the sum of the terms | Y^i P =i 9^ [9ji,p + 9jpj) > does n °t transform as an affine con- 
nection under the TCT. This condition is, however, satisfied by adding the Levi-Civita term 
— 1/2 Xl^=i # y 5W- Now, if one wants the Universal Criterion of Evolution satisfied also 
when the system relaxes along a shortest path, without imposing a priori any restrictions on 
transport coefficients, an extra term to the Levi-Civita affine connection should be added. It 
can be checked that the most general expression for this extra term is 1/(2<t) YTj=i A A4 X J gi p j 
Hence, the affine connection my be written as 

r< 1 V" \nv ( d9lj -i- d9jp d9lp \] 1 STr^Yi dgip (7\ \ 

^ = 2^i 5 [dX- p + dX-dX^)\ + ^^ XX dX (71) 

By performing the same calculations as in Section 1111} we can easily check that the evolution 
equation (parametrized with g), obtained by taking the covariant derivative of both sides of 

5q. ( I48p with affine connection (1711) . is covariant under TCT and satisfies the UCE. In ref. 

l| it is shown that in the general case, (i.e. when = g^ + /^), the expression of the affine 
connection is given by Eq. (lAljl . which generalizes Eq. (17X1) . The covariant derivative, made 



with the affine connection Eq. (lAlj) [or, wit 
as the thermodynamic covariant derivative 



i Eq. (171 P when Ty is symmetric], is referred to 
II- 
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In the examples examined in this paper, the small parameter e measures the distance of 
the system from the steady state. In the article cited in the footnote [19j, it is analyzed 
the relaxation of magnetically confined plasmas in a toroidal geometry. The characteristic 
feature of these equations is the presence of a small parameter, the drift parameter, defined as 
the Larmor radius over a macroscopic length 11] . In this situation, the long-time behaviour 
of the solution, describing the evolution of the system near the steady state, may be obtained 
by using the multiple time-scale perturbation expansion (see, for example, the book cited 



in Ref. 



121] ). The results obtained are valid generally because the dynamics include all 



relevant moment equations and the parameter e is not related to the distance of the system 
from the stationary states. It is possible to show that, even in this case, the results (i) and 
(ii) (see the Introduction) remain valid. 

A non-Riemannian geometry may successively be constructed out of the components of the 
affine connections [lj, jl^ . 

The main conclusion of the present analysis is thus: 

Close to steady- states, the geometry of the thermodynamic space is non-Riemannian with 
affine connection given by Eq. ( L4 The knowledge of the expression for the affine connec- 
tion is a fundamental prerequisite for the construction of the (nonlinear) closure theory on 
transport processes. In Ref. [l|, the curvature tensor and the nonlinear transport equations 
have successively been derived from Eq. (lAlj) and by introducing the following assumption: 
There exists a thermodynamic action, scalar under TCT, which is stationary with respect 
to arbitrary variations in the transport coefficients and the affine connection. 
From this principle, a set of closure equations, constraints, and boundary conditions have 
been derived. These equations determine the nonlinear corrections to the linear (" Onsager" ) 
transport coefficients. The validity of the nonlinear transport equations has been largely 
tested by analyzing several symmetric transport processes, such as the thermoelectric effect 
and magnetically confined plasmas jl4| . Transport processes in tokamak plasmas are of 
particular interest. In this case, even in absence of turbulence, the state of the plasma 
remains close to (but, it is not in) a state of local equilibrium. This prevents the transport 
relations from being linear. The nonlinear transport equations have been used for computing 
the particle and energy losses in tokamak-plasmas and for deriving the distribution density 
functions for the species a in several collisional regimes 
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Appendix A: The Expression of the AfRne Connection and its Main Properties. 



In this section we shall prove the validity of Eqs ( 153]) and Eq. fl55|) . The expression of the 
affine connection is reported in Ref. [l] 



p,q=l 



Q 



I Pi<7i s = 1 



^ L 2(7 yp dXi 2a yp V 



dXi' dXi 



+ ^ J ps A A { dxq , + 

where the thermodynamic Christoffel symbols of the second kind are introduced 

dg jq/ 



{jq j P '=i 



qp > ( dg P 'j dg plq 

dXi' dXi dXP' 



p,g=l 



9ij — 2 ( r «i T i i ) 



P =i 



fij — 2 ( r *J T W 



a denotes the entropy production of the system 

j r T 

a = ^J2 T ^ A o = E T ^ xixi 

i,j=l i,j=l 

We can check that matrices iV u , gij and /y possess the following properties 

r r r r 

a=J2TijX i X i =J29ii Xixj ; E ./o-V'.V ; » ; A"' .V 

ij'=l ij=l i,i=l P=l 



(Al) 



(A2) 



and matrices iV u , #y and are defined as 

r r 

N ij = g ij + - y £ / (f ip g jq + f jp g iq )X*X« with A iZ : ^ .Y ; '*.Y,.,, 8) (A3) 



(A4) 



(A5) 
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Multiplying Eq. (lAlj) with the factor Jj dX^ jdg dX p ' jdg and summing over the indexes, we 
obtain Eq. (|55|) 

^ .dPffA dXidXP ^ dXUXP ^ „„ dXidX p 



i.p,<?=i 



j,p,<?=i 



(A6) 

where identities ( 1A5I) have been taken into account. 

Let us now prove that h % ~ 0(e) and ^I=i ~ ^( e2 ) where e = (f? — l)/l -C 1. The 
definition of /i l is reported in Eq. ( 152]) 



i.p.9 =1 i,p,p' i 9=i 
Taking into account definition (1A3|) and expansions fj4Tl) . we have 

E; =1 fe^f = ^ = 2<#>(e - /) + 42^> - Z) 3 + 0(e») 

I E; =1 iW? = f = 2 4?V - + 42 4^ - S + O(e^) 
Considering that dX % jdg expands as 

dX i 



Xt 1)+ 3XUg-iy + 0(e*) 



dg -d)-— (3JV 
the terms of the second expansion in Eqs. (1A8I) are balanced if we have 

dX i r/.Y-'' 



9ij,p~ 



Hence, close to the steady-state, we have 

r 

a = a st . + 0(e) ; a st . = J] g^X^X, 
1 



JO), 1 Cf(0)_(0) , f(0)_(0Wp y 

% + 7T~ 2^ UiP 9 Jl + Jjp 9iq > X st X , 
P,9=l 



+ 0(6) 



By inserting expansions ( l4Tj) . (1A8|) . (1A10|) and (1A11[) into Eq. f lAlj) . we obtain 

r }p = r }p(l)(^- / )+0(6 3 ) 

where r*. p ^ is the expression flAlj) evaluated at the lowest order. Noting that 

r 

(Q^M^Jjlp g" x M™J q = -6Xi 1} (g - I) + 0(e 3 ) 

j,p,q=i 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 



20 



From expansions (jHJ), (1A12j) and (1A13I) . we finally get 



h l = - 
Thus, at the lowest order 



6X h + E r ki) X (i) X fi) 
J')P=1 



^-/) + 0(e 3 ) 



~ 0(e) 

Multiplying Eq. (|52[) by Jj and summing over the index "i", we obtain 



i=l 



(A14) 



(A15) 



(A16) 



where Eq. (1A6|) has been taken into account. Eq. (IA16I) may also be rewritten as 

r , r r 

E J ^ = 'Jo [E J ^ 1)X W^ " + + E ^(^.)^i)^i) + 0(e 2 ) (A17) 

2=1 



1=1 



E ^ 1)x (d + E ^(^o^D^D + o( e 2 



i=i 



However, from Eq. ( 131]) and Eq. ( |4~3]) . and taking into account that 2~^I=i ^iPCsO^Qi) — 0, 
we have 



E J * 1)x h = E ^(^,)^d4) + o(6 2 ) 

i=l m=i 
Hence, at the lowest order, we finally get 



(A18) 



(A19) 



i=i 
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